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ABSTRACT 

Aims. Gas outflows are believed to play a pivotal role in shaping galaxies, as they regulate both star formation and black hole 
growth. Despite their ubiquitous presenee, the origin and the acceleration mechanism of such powerful and extended winds is not yet 
understood. Direct observations of the cold gas component in objects with detected outflows at other wavelengths are needed to assess 
the impact of the outflow on the host galaxy interstellar medium (ISM). 

Methods. We observed with the Plateau de Bure Interferometer an obscured quasar at z~1.5, XID2028, for which the presence of an 
ionised outflow has been unambiguously signalled by NIR spectroscopy. The detection of ^^CO(3-2) emission in this source allows 
us to infer the molecular gas content and eompare it to the ISM mass derived from the dust emission. We then analyze the results in 
the context of recent insights on scaling relations, which describe the gas eontent of the overall population of star-forming galaxies at 
a similar redshifts. 

Results. The Star formation efficiency (~ 100) and gas mass (Mgas = 2.1 - 9.5 x 10^^ M©) inferred from the CO(3-2) line depend on 
the underlying assumptions on the excitation of the transition and the CO-to-H2 conversion factor. However, the combination of this 
information and the ISM mass estimated from the dust mass suggests that the ISM/gas content of XID2028 is significantly lower than 
expected for its observed M*, sSFR and redshift, based on the most up-to-date ealibrations (with gas fraction <20% and depletion 
time scale <340 Myr). 

Conclusions. Overall, the constraints we obtain from the far infrared and millimeter data suggest that we are observing QSO feedback 
able to remove the gas from the host. 
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1. Introduction 

There are both theoretical (e.g. Hopkins et al. 2008) and obser¬ 
vational (e.g. Sanders et al. 1988; Yan et al. 2010) arguments that 
support the notion that luminous star-forming galaxies (here¬ 
after: ‘Starbursts’) and luminous, unobscured Active Galactic 
Nuclei (AGN; hereafter luminous AGN or ‘QSO’) are basi¬ 
cally the same systems caught in different stages of the co-eval 
growth of (massive) galaxies and the Super Massive Black Holes 
(SMBH) sitting in their centres. In particular, Starbursts should 


Based on observations with the Plateau de Bure millimetre interfer¬ 
ometer, operated by the Institute for Radio Astronomy in the Millime¬ 
tre Range (IRAM), which is funded by a partnership of INSU/CNRS 
(France), MPG (Germany) and IGN (Spain). 
email:marcella.brusa3 @unibo.it 


trace objects caught in the rapid SMBH growth phase charac¬ 
terized by efficient Star Formation (SF), in a dust-enshrouded, 
dense environment, while the unobscured QSOs are systems ra¬ 
diating at the Eddington limit, where the SMBH is almost fully 
assembled. 

Given that both SF and AGN activity are thought to be sus¬ 
tained by the availability of cold gas in galaxies (see e.g Menci 
et al. 2008; Vito et al. 2014), millimeter observations of molecu¬ 
lar transitions are needed to directly probe the presence and state 
of this gas. In the past decade, observations of cold molecular 
gas reservoirs at high redshift (see Carilli & Walter 2013 for a 
comprehensive review) turned out to be crucial in studying the 
gas content and consumption rate in both normal and peculiar 
systems. For example, the gas properties of “normal" galaxies 
are being investigated in increasing details up to high-z (Tac- 
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coni et al. 2013; Genzel et al. 2015; Sargent et al. 2014), and as 
a function of many of the structural and physical properties of 
the systems (e.g. Star Formation Rate, SFR; stellar mass; colors; 
see e.g. Genzel et al. 2015; Sargent et al. 2014). This has be¬ 
come possible thanks to the large investment of time at millime¬ 
ter arrays, mainly the Plateau de Bure Interferometer (PdBI). 
In particular, it has been reported that, among massive systems, 
(M^ > 10^^ Mq), the gas fraction increases across the main se¬ 
quence (MS; defined between the SFR and the stellar mass of 
galaxies) at fixed redshift (see Magdis et al. 2012a,b; Saintonge 
et al. 2012; Tacconi et al. 2013; Sargent et al. 2014) and is hence 
closely related to the Specific Star Formation rate (sSFR). This 
fits in a scenario where the redshift evolution of the sSFR is 
consistent with being driven by the gas fraction (see also Lilly 
et al. 2013). Similar conclusions are reached in works involving 
dust fitting methods to derive the gas mass (see e.g. Santini et al. 
2014). 

The first molecular studies on local Ultra Luminous In¬ 
frared Galaxies (ULIRGs) and Submillimeter Galaxies (SMG) 
at higher redshifts, i.e. targetting objects in the ‘Starburst’ phase, 
showed that these systems typically have a low molecular gas 
content with respect to their current SFR, or alternatively higher 
star-formation efficiencies. Indeed, defining empirically the Star 
Formation Efficiency (SEE) as the ratio of the IR luminosity 
to the CO luminosity (in units of Lo/(K km s“^ pc^)), ‘Star- 
bursts’ have SFE>200 (see e.g. Daddi et al. 2010; Genzel et al. 
2010) larger than those observed in normal star forming galax¬ 
ies with the same molecular gas content (Tacconi et al. 2010, 
SEE~ 50 - 200). In other words, their consumption time scale 
is shorter with respect to normal galaxies and they will exhaust 
their gas reservoirs in a short timescale (<100 Myr). This is 
consistent with the hypothesis that ‘Starbursts’ in general (and 
ULIRGs/SMGs in particular) are objects at the peak of their SE 
activity in the heavily obscured phase. 

On the other hand, high values of the Lir/L’ (CO) ratio have 
been also observed in high-z unobscured QSO host galaxies 
(SEE> 200; e.g. Solomon & Vanden Bout 2005; Riechers 2011; 
Riechers et al. 2011), although, being a subsequent phase of 
‘Starbursts’ in the evolutionary sequence, their SER is expected 
to be already substantially suppressed. In this case a significant 
fraction of the gas could have been previously removed during 
the ‘blow-out’ phase, and the observed high SEE in unobscured 
QSOs can be ascribed to region of residual, on-going SE, point¬ 
ing towards a possible effect of ’positive feedback’ on the galaxy 
from the AGN (Silk 2013; Zubovas & King 2014). 

What is still missing for a full understanding of the results 
of the aforementioned studies, in terms of the role of the phys¬ 
ical processes which govern the co-eval BH-galaxy growth, is 
a full characterization of the gas properties of objects caught 
in the short-lived “transition” phase between the Starburst and 
QSO stages. This phase is expected to be characterized by gas 
reservoirs not yet depleted and by complex kinematics, includ¬ 
ing strong winds and outflows. Brusa et al. (2010) proposed that 
sources in the ‘blow-out’ phase at z~ 1.5 can be isolated on the 
basis of their observed X-ray-to-optical-to-NIR colors and pre¬ 
sented the source XID2028 (z= 1.5927), detected in the XMM- 
COSMOS survey, as the prototype of this class. XID2028 is a 
luminous (Lboi -2x10^^ erg s“Q, mildly obscured QSO hosted 
in a massive galaxy, with M* ~ 4.5 x 10^^ M© and a SER~ 270 
M© yr“^ as measured by Herschel from PEP and SPIRE data 
(Lutz et al. 2011; Bethermin et al. 2012). At its center, XID2028 
has a supermassive black hole with mass Mbh ~ 3 X 10® Mo 
(Bongiomo et al. 2014), which is accreting at ~ 5% of its Ed- 
dington luminosity. 


The presence of a massive outflow in the ionized gas com¬ 
ponent of XID2028, traced by the [O III]/15007 emission, has 
been unambiguosly and independently confirmed by X-shooter 
slit spectroscopy (Brusa et al. 2015; Perna et al. 2015) and SIN- 
EONI J-band IFU observations: in fact, XID2028 hosts one of 
the most massive {Mion > 250 M© yr“\ with v> 1500 km s“Q 
and most extended (out to scales of ~ 13 kpc) outflows detected 
in a high-z QSO (Cresci et al. 2015). Most importantly, the out¬ 
flow lies exactly in the center of a cavity in star forming regions 
in the host galaxy (as traced by narrow Ha emission line map 
and rest frame U band imaging; see Cresci et al. 2015) thus sug¬ 
gesting that the wind is removing the gas from the host galaxy 
^negative feedback’), and at the same time is also triggering star 
formation by outflow induced pressure at the edges (‘positive 
feedback’; e.g. Zubovas & King 2014). XID2028 therefore rep¬ 
resents a test case to study QSO ‘feedback in action’. However, 
the evidence of feedback in this source mostly comes from mea¬ 
surements of the on-going star formation in the source traced by 
the narrow Ha emission line that in principle may be affected by 
e.g. differential extinction effects in the host galaxy. 

Direct observations of the cold gas component in this galaxy 
are needed to assess whether the ionized outflow has an impact 
on the cold gas reservoir. With this aim, here we present ob¬ 
servations of the CO(3-2) transition of XID2028, redshifted to 
2mm, obtained with the PdBI Interferometer. We compare the 
gas masses derived from CO with that inferred from the dust 
mass and based on Far Infrared (FIR) data. These two meth¬ 
ods allow us to investigate whether AGN feedback has already 
been effective in diminishing the cold gas mass in the host, or 
whether the feedback phase is still associated with cold-gas-rich 
galaxies similarly to MS star-forming galaxies, with important 
consequences for galaxy-AGN coevolutionary models. 

The paper is organised as follows: Section 2 presents the 
PdBI observations and data analysis. Section 3 discusses the re¬ 
sults, while Section 4 summarizes our conclusions. Throughout 
the paper, we adopt the cosmological parameters //q = 70 km s“^ 
Mpc“\ Q;„=0.3 and Qa= 0.7 (Spergel 2003). In quoting magni¬ 
tudes, the AB system will be used, unless otherwise stated. We 
adopt a Chabrier Initial Mass Function to derive stellar masses 
and SFRs for the target and comparison samples. The physical 
scale is F'- 8.5 kpc at the redshift of the source. 

2. Millimeter observations 

2.1. Data reduction 

XID2028 was observed with receivers tuned to a frequency of 
133.37 GHz, corresponding to the expected frequency of the 
CO(3-2) emission line, with the PdBI array in the most com¬ 
pact (D) configuration. The observations were split in 3 tracks 
(31-May, 1, 6, June 2014). The system temperature (T^ 3 ;^) was 
between 100 and 300 K, and water vapor 4-6 mm. The quasar 
1005-r058 (0.3 Jy at 133.7 GHz) was used as a phase and am¬ 
plitude calibrator. MCW349 (with a flux of 1.8 Jy) was used for 
absolute flux calibration, which yields an absolute flux accuracy 
of about 5% at the observed frequency. Calibration and map¬ 
ping were done in the GILDAS environment. The flagging of 
the phase visibilities was fixed at < 45° rms. 

The total observing time was 5.6 hrs (3.06 hrs on source), for 
a total of 3673 visibilities available, before applying any flag. We 
then removed one scan (3994 in 01-June track) due to problems 
with the tracking. Data from antenna 1 from the 06-June track 
were not used in the final dataset due to the presence of a tuning 
parasite that produced a spurious signal at a frequency (133.21 
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Fig. 1: Left panel: spectrum of XID2028 integrated over the beam. The solid line shows a Gaussian fit with FWZI=770 km s“^ 
and centered at the frequency corresponding to the redshift of the source. Right panel: integrated map of CO(3-2), in the channels 
corresponding to the “systemic" peak of the line. Contour levels are Icr each (cr=0.23 Jy km s“^). The synthesised beam is shown 
in the bottom-left corner. The black cross marks the phase center (i.e. the ACS nucleus). The blue and red cross mark the positions 
of the blue and red line components, as derived from our spectroastrometric analysis. 


GHz) close to the observed frame frequency of the CO(3-2) tran¬ 
sition. After flagging of bad visibilities, the total on source time 
is 2.54 hours (six-antenna equivalent), and the Icr sensitivity is 
1.36 mJy/beam in 20 MHz channels, for a total of 3052 visi¬ 
bilities. The clean beam of the observations is 4.5"x3.4", with 
an angle of 38 degrees. The phase center of the data set was 
set to the HST position of the QSO nucleus (RA=10:02:11.29, 
DEC=-r01:37:06.79). 

2.2. Analysis 

We estimated the 2 mm continuum by collapsing the line-free 
channels of the data set and fitting the visibilities. The continuum 
is not detected with a 3cr upper limit on its flux of 0.3 mJy. 

The redshift of the host galaxy (z= 1.5927) was adopted to 
convert the frequency to velocity space. Figure 1 shows the line 
spectrum integrated over the beam. The line displays two peaks: 
one centered around the systemic redshift (FWHM~550 + 200 
km s“^ from a Gaussian fit), and another centered at ~ 1000 km 
s“^ (henceforth referred to as “red feature"). The peak at the sys¬ 
temic position is significant at 5cr, while the “red feature" is at 
a lower significance (~3 (t). Moreover, the “red feature" peaks at 
~ 133.0 GHz, close to a known parasite signal at 132.9 coming 
from antenna 4 and identified in all tracks. We created a new 
table flagging data on Antenna 4. The total exposure time de¬ 
creased to 1.8hr and the total number of visibilities to be used 
for the scientific analysis also considerably decreased. The red 
feature is not significant anymore (S/N< 3). However, the sig¬ 
nificance of the detection over the systemic line also decreased 
at S/N~4. For this reason we decided to keep the full datasets in 
the analysis and, in the absence of deeper, more highly resolving 
observations which could confirm the presence of a second dy¬ 
namically distinct component, we will consider the red feature 
as spurious. 

The zero spacing flux estimated by fitting the averaged vis¬ 
ibilities in the velocities range from -340 to -1-440 km s“^ with 
a point source function is S’(CO)=1.6±0.3 mJy (5.3cr), and re¬ 
turns a centroid at (RA,DEC=10:02:11.24 01:37:05.48). The in¬ 
tegrated flux over the full velocity range of the systemic line 


(with Full Width Zero Intensity, FWZI~ 770 km s“^) is there¬ 
fore f S'(CO)dv= 1.23+0.23 Jy km s“^ This measure depends 
only on the data calibration (including flagging of the Antennas) 
and does not depend on any other assumption, like e.g. mask¬ 
ing, extraction region, ad-hoc centroid. The quoted errors take 
into account the statistical errors of the nv plane fit and the er¬ 
rors on the absolute flux calibration (5%). The right panel of 
Figure 1 shows the integrated map over the systemic line emis¬ 
sion. We verified that the flux extracted from the integrated map 
(8=1.55+0.3 mJy) on a region slightly larger than the beam, is 
in agreement with the one estimated by fitting the visibilities. 

Figure 2 shows the HST/ACS image (background) with su¬ 
perimposed the contours from the K-band image (blue), which 
should trace the extension of the host galaxy. The black contours 
are from the map obtained on the line detected at the systemic 
position (e.g. from Figure 1 right, in steps of S/N, starting from 
Icr) and the black cross marks the line centroid. 

From both Figures 1 and 2 it is clear that the line peak is 
offset by ~ 1" from the QSO nucleus position. From previous 
observations with the same phase calibrator (1005+058), we can 
exclude errors in the absolute astrometry. The error associated 
with the beam and the S/N of the source translates into a po¬ 
sitional uncertainty of 0.46" x 0.36". We note, however, that 
the displacement may be due to the limited nv coverage of the 
data, and that a CO-offset is typical of low S/N data (see e.g. 
Casey et al. 2011). Better signal-to-noise ratio and nv coverage 
are needed to refine the location of the gas reservoir. 

A dynamical mass can be estimated from the CO line width 
assuming a size (R) and an inclination (/) of a rotating molecular 
gas disc. The size can be inferred using the spectroastrometric 
technique (Gnerucci et al. 2011; Camiani et al. 2013, and refer¬ 
ences therein), applied to the CO data cube. By integrating the 
CO data in the red (0,+400 km s“Q and blue (-400,0 km s“Q 
line channels, we measure a difference in the line centroids of 
~ 1.5 + 0.2" (with an error of 0.3 pixels for each detection). The 
centroids of these detections are also shown in Figure 1 (right 
panel) as blue and red crosses to mark the blue and red line chan¬ 
nels, respectively. The measured shift corresponds to ~13 kpc at 
the source redshift and translates to R~6.5 + 0.8 kpc, in agree- 
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Fig. 2: HST/ACS image (F814W filter) with superimposed K- 
band contours from CFHT (blue, arbitrary levels chosen to trace 
the whole K-band emission). Black contours represent CO(3-2) 
emission from the integrated map in the channels corresponding 
to the “systemic” peak of the line (same levels as in right panel 
of Figure 1; starting from 1 sigma). The black cross (with asso¬ 
ciated ellipse) marks the line centroid. The image is about 10" 
across. The beam size is 4.5"x3.4", with an angle of 38 degrees. 

ment with the extension seen in the K-band data (see also Cresci 
et al. 2015). Applying Equation 5 of Gnerucci et al. (2011), we 
infer a Mdyn(sin i)^=4.5xl0^^ M© and, assuming an inclination 
of 60 deg, a Mdyn~6.0+2.3 x 10^^ M© once all the uncertain¬ 
ties in the quantities are taken in account. We will discuss in the 
following how this compares with the total mass derived from 

M-*, -|- Mgas • 

3. Results and discussion 

Deriving the luminosity L’CO [1-0] of the ground-state transition 
(which is generally regarded as the best indicator of the total gas 
reservoir) requires an assumption on r 3 i, the luminosity ratio be¬ 
tween the CO(3-2) and the CO(l-O) transitions, which depends 
on the nature of the systems: a ratio of ~ 0.7 - 1 is typically re¬ 
ported for SMG galaxies and QSOs (see Carilli & Walter 2013 
and references therein), while an average ratio of ~ 0.42 has been 
determined for MS star-forming galaxies at a similar redshift as 
XID2028 (e.g. Daddi et al. 2015). 

The CO(3-2) luminosity of XID028 is L’(CO[3- 
2])=1.9+0.4 X 10^^ K km s”! pc^ (~ 9.3 + 2 x 10^ L©), 
following Solomon & Vanden Bout (2005). This value is in 
between the average values observed in this molecular transition 
for U/LIRGs (L[CO(3-2)]=2.6 + 0.5x10^ K km s'^ pc^) and 
SMGs (L[CO(3-2)]4.4 + l.lxlO^^ K km s-^ pc^), as reported 
in the work of lono et al. (2009). On the other hand, the SFR 
(~ 270 M© yr-i) and M* (~ 4.5 x lO^^ M©) of XID2028 are 
consistent with those observed in a MS galaxy at z~ 1.5 (see 
Mainieri et al. 2011; Bongiomo et al. 2014; Brusa et al. 2015). 

The CO(2-l) transition in XID2028 is not detected down 
to a sensitivity of 0.23 mJy/beam over the full 770 km/s line 
width (corresponding to a 3(T upper limit on the line integrated 
flux of 0.53 Jy km/s), from a separate, 3mm-band PdBI ob¬ 
servation of XID2028 in October 2014 (M. Sargent, private 


communication). This suggests a near thermal CO-excitation 
stated and therefore rsi around unity, i.e. larger than the standard 
value usually adopted for MS galaxies, more consistent with the 
QSO/Starburst scenario. Given the complex nature of the sys¬ 
tem, we derive the CO(l-O) luminosity under the conservative 
assumption that r3i=0.7 (consistent with the constraints we have 
from millimeter data alone), and we will apply different aco 
factors to derive molecular gas masses under the QSO/ULIRG 
and MS assumptions, discussing the implications of the findings 
in the different cases. The inferred L’(CO[1-0]) luminosity for 
XID2028 (abbreviated as L’(CO) in the following) is therefore 
L’(CO)=2.6xlOi^ Kkm s-i pc^ (~ 1.2 x 10^ L©). 


3.1. Star Formation Efficiency 

Figure 3 (left panel) shows L’(CO) against the total Infrared Lu¬ 
minosity (Lir, computed between 8-1000 pm) for XID2028 (red 
circle). The IR luminosity of XID2028 is very well constrained 
by Herschel/PACS-rSPIRE data (logLiR=12.47; see Brusa et al. 
2015; Perna et al. 2015) and has been estimated from fitting 
all bands with photometry with rest-frame wavelength > 50yum 
with Dale & Helou (2002) Starbursts templates, using the same 
technique as in Santini et al. (2009). Although recent works on 
FIR emission of AGN show that even the flux observed at rest 
frame wavelengths longer than 60yum can be AGN-dominated 
(e.g. Mullaney et al. 2011), in XID2028 the QSO contribution 
is expected to be negligible, as shown in Figure 2 of Pema et al. 
(2015), where the most recent SED fitting decomposition for this 
object is presented. The observed IR luminosity corresponds to 
a SFR of ~ 270^ jQQ M© yr“^ using the SFR-IR luminosity rela¬ 
tion (Kennicutt 1998), and taking into account the uncertainties 
on the flux normalization of the starburst component related to 
the AGN-host SED decomposition. 

We compare this measurement with the compilation of low- 
and high-redshift normal star-forming galaxies with measured 
aco presented in Sargent et al. (2014) and the SMG-sample from 
Bothwell et al. (2013), involving both outliers and galaxies con¬ 
sistent with the locus of the main sequence at their redshift. We 
also plot the unobscured QSOs at 1< z <4 presented in Riech- 
ers (2011)^. For these sources, the IR luminosities are extracted 
from the Carilli & Walter (2013) compilation. Finally, in Fig¬ 
ure 3 we show SW022550 and SW022513 at z~ 3.4 (Polletta 
et al. 2011), ULASJ1539 at z~ 2.5 (Feruglio et al. 2014), and 
the MIPS-selected sources at z~2 from Yan et al. (2010). All 
these systems have been proposed to be in the “transition phase" 
between an heavily obscured Starburst phase and the unobscured 
QSO phase. 

The SEE of XID2028 (SEE- 110) is on the lower side 
of the SFEs measured for high-z SMG and unobscured QSOs 
(SEE- 100 - 1000). Instead, the SEE is consistent with those 
reported (albeit with much larger uncertanties due to the lack 
of a complete multiwavelength coverage and reliable measure¬ 
ments of Lir) in the obscured QSOs systems proposed to be in 
the “transition phase" mentioned above. 


^ r3i>1.0 assuming CO[2-l] is thermalized, and r3i20.9 assuming 
CO [2-1] is sub-thermally excited r 2 i = 0.84, standard value for MS 
objects 

^ In the case of lensed quasars, the values are corrected for the ampli¬ 
fication, as reported in Riechers 2011. 
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Fig. 3: Left panel. L’(CO[1-0]) against the IR luminosity (S-lOOOyum) showing a compilation of MS galaxies at 0<z<2.5 from the 
Sargent et al. (2014) work (Sal4; grey symbols), local ULIRGs from from , So97, crosses SMGs from Bothwell et al. (2013; Bol3, 
filled brown squares), l<z<4 QSOs with data from Riechers (2011) and Carilli & Walter (2013; R11/C&W13, light blue circles). 
The measurement for XID2028 obtained assuming rai = 0.7 is shown as a red circle. Obscured QSOs proposed to be in the transition 
phase presented in Aravena et al. (2008), Polletta et al. (2011), Feruglio et al. (2014) and the MIPS-selected sources from Yan et 
al. (2010) are also marked, as labeled. Right panel Inverse, integrated Kennicutt-Schmidt relation between SFR and molecular gas 
mass. The color points correspond to the position of XID2028 with different CO-to-H 2 conversion factors: afco=3.6 (top; blue), and 
aco=0.8 (bottom; red). The green star shows the ISM mass inferred from the dust SED. All the values for XID2028 are slightly 
offset in the x-axis for clarity. Other points are taken from the same samples presented in the left panel. For the Bol3, R11/C&W13 
and YalO samples we plot the median value with associated 16% and 84% percentiles. For the obscured QSOs all the authors used 
In both panels, the solid black line is the best-fit relation for MS galaxies, the dashed line defines the locus of strong SB 
galaxies with approximately 15 times shorter depletion time (Mmoi/SFR) than MS galaxies (Sargent et al. 2014). 


3.2. Molecular gas mass from CO data 

Estimating the molecular gas mass based on the CO luminosity 
critically hinges on the CO-to-H2 conversion factor aco, defined 
as the ratio between the mass of molecular gas (Mmoi) to the in¬ 
tegrated CO(l-O) luminosity (aco=Mnioi/L’(CO)[l-0] in units of 
Mo/(K km s“^ pc^)). This value depends on the ISM conditions, 
and two distinct assumptions are often adopted: aco ~ 4 for ex¬ 
tended SF disks/MS galaxies of solar metallicity, and aco ~ 0.8 
for compact luminous systems (Downes & Solomon 1998; see 
Carilli & Walter 2013 and Bolatto et al. 2013 for in-depth dis¬ 
cussions). 

From a morphological point of view we do not have a clear 
classification on the properties of the host galaxy. Given that the 
HST image suffers from substantial extinction (Ay ~ 3; see dis¬ 
cussion in Pema et al. 2015) and, in any case, is dominated by the 
central active nucleus, it cannot be used for a reliable morpho¬ 
logical analysis. However, no clear signatures of merging struc¬ 
tures are visible in the rest-frame U-band. The low-resolution 
(with respect of HST) K-band image is instead consistent with 
both an elliptical galaxy and a spiral galaxy, possibly interacting 
with a north-east system (see Figure 2). Even if the MS is mainly 
populated by "normal" spiral and disk galaxies (see e.g. Wuyts 
et al. 2011a), we note that XID2028 would lie among the popu¬ 
lation which occupies the upper envelope of the MS at z~ 1.5. 
These galaxies may have also cushier light profiles, intermedi¬ 
ate between disky galaxies and red and dead systems (see Wuyts 
et al. 2011b, their Figure 1, right panel). In any case, if after 


point-source subtraction this galaxy were to be shown to have 
an early-type or disturbed host galaxy morphology, it would ac¬ 
tually be highly consistent with the statistical findings of Wuyts 
etal. (2011b). 

In the vast majority of studies targeting SMGs, QSOs and 
ULIRGs systems (see e.g. Aravena et al. 2008; Riechers et al. 
2011; Polletta et al. 2011; Feruglio et al. 2014, among others), 
crco = 0.8 has been adopted even in absence of better informa¬ 
tion on the physical properties of the system (e.g. compactness 
of the source). Under the assumption of starbursts/QSO scenario, 
we obtain for XID2028 a gas mass Mmoi ~2.1+0.4 x 10^^ M©. 

To infer the molecular gas mass under the MS hypothesis, 
we consider a metallicity dependent conversion factor aco (see 
e.g.Genzel et al. 2012; Bolatto et al. 2013). In the following we 
will assume for XID2028 a value of 12 + logiOjH) = 9.07, 
the metallicity inferred from the so-called Fundamental Metal¬ 
licity Relation (FMR, Mannucci et al. 2010), that relates the 
metal content with the stellar mass and the SF of the galaxy inde¬ 
pendently of redshift (Cresci et al. 2012). Applying the relations 
describing redshift-dependent variations of aco in the SFR-M* 
plane of Sargent et al. (2014) to XID2028 one would expect^ 
aco ~ 3.6, and the corresponding molecular gas mass would 
then be M^ioi ~ 9.5 ± 1.9 x 10^^ M©. 


^ A virtually identical conversion factor would be inferred using the 
relation between metallicity and aco calibrated in Genzel et al. (2012) 
once the offsets between different metallicity calibrations are taken into 
account. 
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The two values for Mgas inferred under the two different as¬ 
sumptions are plotted in Figure 3b (with the statistical errors as¬ 
sociated to the line detection), where the molecular gas mass is 
shown as a function of the SFR (for the same samples presented 
in Fig. 3a). 

3.3. Molecular gas mass from FIR emission 

We also adopt an independent method to compute the total gas 
mass in this source, using the dust mass derived from FIR pho¬ 
tometry. For this purpose, we assume a metallicity-dependent 
gas to dust ratio, following the calibration presented by Santini 
et al. (2014) and recently extended to AGN samples in the work 
by Vito et al. (2014). This estimate is independent from aco^ al¬ 
though it depends on the metallicity (Z) of the system and on 
the assumptions that the dust-to-gas ratio scales linearly with Z 
through a constant factor (Draine et al. 2007). 

The dust mass is obtained via SED decomposition of the 
AGN and host galaxy contributions, using a combination of the 
Silva et al. (2004) AGN templates and the Draine & Li (2007) 
dust templates to fit the 100-500yum range. For objects at z>l, 
submillimeter data are in principle required to properly sample 
the dust emission free from AGN contamination. However, we 
note that for XID2028 the best fit SED decomposition performed 
following Vito et al. (2014) is consistent with the upper limit of 
the continuum at 2mm (see Section 2.2). 

We obtain a total dust mass of Mdust = 7.7 + 4.2 x 10^ M©. 
Assuming the FMR metallicity (see above) and the Santini et al. 
(2014) calibration, this translates into a Mgas ~ 4.5 + 2.4 x 10^^ 
M© (without considering the uncertainty on the dust-to-gas-ratio 
calibrations, e.g. a factor of ~ 2, see Sandstrom et al. 2013). We 
note that using the Leroy et al. (2011) metallicity dependence 
of the dust-to-gas ratio would yield consistent results, within the 
errors. 

The value inferred from the dust fit approach is plotted as a 
green star in Fig. 3b. If we use this estimate for Mgas, and the 
observed L’(CO), we can derive an effective aco for this source, 
trco(dust) ~ 2.4 X rsi (aco(dust) ~1.7 given our adopted excitation 
correction r3i=0.7). 

3.4. Gas fraction and depletion timescale 

The uncertainty in the derived gas mass from the CO data is 
dominated by the assumption in aco (a factor of 4.5) with re¬ 
spect to the statistical uncertainties (20%). Given that the value 
derived from the dust fit is in between those for the two different 
assumptions in aco^ m the following we will refer to this value as 
our best estimate for the molecular gas mass, and those from the 
CO data under the MS and QSO/Starburst assumptions as upper 
and lower limit, respectively, i.e. Mgas = 4.5(1.7 - 11.4) x 10^^ 
M©, where the lower and upper limits in parenthesis take also 
into account the statistical uncertainties of the detection, and 
overall also the uncertainty in the assumed dust-to-gas ratio. We 
note that the total gas mass inferred from the dust continuum fit 
includes both the molecular and atomic components. However, 
the atomic mass usually constitutes a negligible fraction of the 
total gas mass. 

The stellar mass of XID2028 is M^ ~ 4.5 x 10^^ M© from the 
most recent SED fitting decomposition (Pema et al. 2015). This 
value is a result of the inclusion in the multicomponent SED fit¬ 
ting of a mildly obscured QSO component, given that we observe 
the Broad Line Region (BLR) emission in the Ra line complex 
(Bongiomo et al. 2014). In Section 2.2 we reported a dynamical 


mass Mdyn~6+2xl0^^ M©. Although the estimate of the dynami¬ 
cal mass suffer from large uncertainties, it is quite reassuring that 
it is consistent with the value we obtain from the sum of the stel¬ 
lar and molecular mass components (Mtot~ 4.6 - 5.6 x 10^^ M©, 
taking into account the range of Mg^^). 

We can then calculate the molecular gas fraction, yUmoh de¬ 
fined as the ratio of the molecular gas mass and the stellar mass 
Oumoi=Mnioi/M*; see e.g. Sargent et al. 2014; Genzel et al. 2015). 
Given the molecular gas masses inferred in the previous Section, 
the gas fraction translates into ~ 5% for the QSO/Starburst and 
~ 21% for the MS scenarios. The value from the dust mass mea¬ 
surement is in between these two estimates (~ 10%). Similarly, 
we can estimate the depletion time scale (defined as Mgas/SFR; 
e.g. the rate at which the gas is converted into stars) and we in¬ 
fer tdepi=75, 340 and 160 Myr using the QSO/Starburst, MS and 
dust-fit derived gas masses, respectively. 

3.5. Evidence for QSO feedback 

Figure 4 (left panel) shows the gas fraction in XID2028 
for the three assumptions described above, plotted against 
the sSFR-excess with respect to the main sequence, e.g. 
sSFR/sSFRms=0-86, where the mass- and redshift-dependence 
of the characteristic sSFR of MS galaxies follows the calibration 
in (Sargent et al. 2014) which is based on a large compilation 
of literature data. In this plot we show the same samples used in 
Figure 3 (with the exception of unobscured QSOs and the MIPS 
selected sources for which no stellar mass estimates are avail¬ 
able) and we plot as a solid line the median trend with normal¬ 
ized sSFR, expected for a galaxy of the same mass and redshift 
of XID2028 (taken from the 2-Star Formation Mode description 
of normal and starbursting (off-MS) galaxies in Sargent et al. 
2014). 

Taking the best Mgas estimate for our target, and even tak¬ 
ing into account the uncertainty on aco assumption, XID2028 
is among the objects with the lowest gas fraction for its sSFR 
detected so far in the high-z Universe and associated to normal 
star-forming galaxies (green star in Figure 4), especially when 
compared to systems with similar masses (solid line). The yUmoi 
is instead more similar to that expected for ‘Starburst’ galax¬ 
ies of a similar mass and redshift (see value of black trend line 
at sSFR/<sSFR>ms^4), but XID2028 does not share with these 
sources the same burst of star formation. Instead, the gas fraction 
of XID2028 is similar to normal galaxies in the local universe 
(open triangles), despite its higher redshift. 

An alternative way of visualizing the gas content and con¬ 
sumption is illustrated in the right panel of Figure 4, where the 
depletion time scale is plotted against the MS-normalised sSFR 
of the host galaxy. Assuming our best Mgas estimate, XID2028 
lies at shorter depletion time scales with respect to MS galaxies 
(at any redshift), i.e. it is consuming its residual gas more rapidly 
than normal star-forming galaxies. This qualifies XID2028 as a 
clear outlier with respect to the average population, and a rare 
object, consistent with the hypothesis that it is caught in the very 
short transition phase in which the QSO feedback is released. 

Similar conclusions can be reached examining the position 
of our source with respect to the Kennicutt-Schmidt relation 
(Kennicutt 1998): assuming the physical scales inferred in Sec¬ 
tion 2.2, and that the molecular gas and the SF episodes are dis¬ 
tributed uniformly over this region, XID2028 would lie slightly 
above (a factor ~ 2.5) the correlation observed for normal and 
starburst galaxies. However, the SFR density measured in this 
way is to be considered a lower limit, given that the SF regions 
seem to be patchy (see Cresci et al. 2015). Therefore, XID2028 
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Fig. 4: Gas fraction yUmoi panel) and depletion timescale (Rightpanel) plotted versus the sSFR-excess for the same samples and 
with the same color code presented in Figure 3. The values for XID2028 are slightly offset in the x-axis for clarity. All quantities are 
normalized to the expected values for normal and starburst galaxies predicted by the calibration presented in Sargent et al. (2014). 
The black line traces the expected variation (median) with sSFR for a MS spiral galaxy with identical mass and redshift as XID2028 
(see Sargent et al. 2014). The step at sSFR/<sSFR>~4 reflects the transition from the main sequence locus to the sSFR-regime 
where high-SFE starbursts dominate. XID2028 lies a factor ~2 to ~10 below the black line, i.e. it shows significant lower gas 
fraction and depletion time scale than those expected for the properties of the its host galaxy. 


would further deviate above the K-S relation, towards regions of 
short depletion timescales. 

It is important to note that, even when using the MS assump¬ 
tion, the molecular gas fraction and depletion timescale would be 
considerably lower than those expected for systems of the same 
host galaxies properties of XID2028 (blue point in Figure 4). In 
particular, the depletion timescale observed in XID2028 for the 
MS scenario is a factor of ~ 2 lower than the expectations of 
Sargent et al. (2014) and a factor ~ 3 lower than that obtained 
by the parameterization of MS and off-MS galaxies presented in 
Genzel et al. (2015, using their global fit we expect for XID2028 
tdep(G15)giobai - 970 Myr). The discrepancy with the calibra¬ 
tions is more extreme if the values obtained in the QSO scenario 
are adopted (red circles in Figure 4). We also note that our chain 
of assumptions in deriving Mgas has been very conservative. For 
example, we used rsi^O.T instead of r3i^.9 as suggested by the 
non detection of CO[2-l] emission, which would have instead 
provided a 20% smaller CO[1-0] flux. This conservative assump¬ 
tion also compensates a possible overestimate of the value of the 
CO[3-2] flux, which could result from measuring the line flux at 
the phase centre rather than at the slightly offset centroid. The 
result of the lack of molecular gas in XID2028, with Mga^lO^^ 
Mo, is therefore quite robust. 

A short depletion time scale with respect to MS galaxies has 
also been found for SMGs in the Bothwell et al. (2013) sam¬ 
ple, and other AGN/Starburst systems plotted in Figure 4 (Ar- 
avena et al. 2008; Polletta et al. 2011; Feruglio et al. 2014). Yan 
et al. (2010) also reported a short depletion timescale of ~40 Myr 
for the sample of MIPS-selected ULIRGs. The short depletion 
time scale in SMGs has been interpreted as higher star forma¬ 
tion efficiency in the galaxy (e.g. Genzel et al. 2010; Daddi et al. 


2010), probably due to higher density of the ISM in these com¬ 
pact systems. These may also be the case for ULASJ1534 and 
the COSBOll, which have sSFRs comparable to SMGs, and for 
which we expect compact gas reservoirs. 

Instead, in XID2028 a significant fraction of the gas is ex¬ 
pected to be already expelled from the galaxy. The SF is then 
probably maintained only in the denser environments, less af¬ 
fected by the negative feedback, and possibly enhanced by pos¬ 
itive feedback due to the outflow induced pressure (e.g. Silk 
2013). The fact that XID2028 has a smaller gas reservoir and 
shorter depletion time than that measured for MS galaxies of 
similar sSFR therefore constitutes a new probe, in addition to 
the analysis presented in Cresci et al. (2015) based on NIR data, 
that QSO feedback in the form of powerful outflows is able to 
affect star formation in the host and expel a significant fraction 
of gas from the host galaxy. 

4, Summary 

We presented the first molecular line luminosity measurement, 
via CO(3-2) observations obtained at the PdBI interferometer, 
in a luminous obscured QSO at z~1.5. The target is thought to 
be in the ‘blow-out’ phase, and the presence of a powerful out¬ 
flow with significant impact on the host galaxy has been unveiled 
through previous NIR observations (Pema et al. 2015; Cresci 
et al. 2015). We complemented the PdBI data with FIR dust fit¬ 
ting, and report the following results: 

• We measure a SFE (^110) at the lower end of those re¬ 
ported in the literature for a large number of QSOs and Star- 
burts/SMG galaxies (see lono et al. 2009; Carilli & Walter 
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2013), and consistent with that inferred for obscured QSOs 
at higher redshift; 

• We infer a molecular gas mass (Mmoi) in the range 2.1 + 0.4 ^ 
9.5 + 1.9 X 10^^ M© applying the QSO/Starburst or MS con¬ 
version factors to the measured L’CO line luminosity, respec¬ 
tively, and a total gas mass Mgas ~ 4.5 x 10^^ M© from dust 
continuum fitting; 

• A value for the molecular gas mass < 10^^ M© is also re¬ 
markably consistent with our estimates of the dynamical 
mass through spectroastrometric methods (see Section 2.2), 
given the high stellar mass of XID2028; 

• We also infer a molecular gas fraction /tmoi ~ 5 - 20%. This 
translates into a gas depletion time scale tdepi -70-340 Myr, 
depending on the assumptions on aco (see Figure 4). 

• The value of tdepi is considerably lower (30%) than those 
observed in systems hosted in similar massive (M^ > 10^^ 
M©) MS galaxies (MS-normalised sSFR~ 1), and consistent 
with those observed for SMGs and for the other few systems 
proposed to be in the transition phase. 

We propose that in XID2028 the QSO wind, detected in the 
ionised gas component out to 10-kpc scales, has already removed 
most of the molecular gas from the host galaxy. All the obser¬ 
vational constraints (low molecular gas content, lowest yUmoi at 
a fixed sSFR when compared to M^ > 10^^ M© systems, and 
lowest sSFR at a fixed Pmoi) are consistent with such a scenario, 
where the gas in the host galaxy of XID2028 is indeed already 
depleted/dispersed by the effects of the strong QSO feedback 
(see also Coppin et al. 2008 and Yan et al. 2010 for similar inter¬ 
pretation). In dense regions (e.g. clumpy M^^^ reservoirs), possi¬ 
bly located at the edge of the outflow cavity (Cresci et al. 2015), 
the residual gas is converted into stars at a high rate similar to 
that observed in SMGs, where the low depletion time scale is in¬ 
deed ascribed to the efficient SF triggered in dense and compact 
gas reservoirs. 

The measure of the intensity of the CO(3-2) emission in 
XID2028 represents a first step towards a mapping experiment 
using high spatial resolution to study the morphology and the 
kinematics of the molecular gas reservoir and of the clumpy 
structures in the distribution of SF regions seen in HST and SIN- 
FONI maps. Sensitive ALMA and/or NOEMA observations of 
XID2028 will finally give the spatial resolution to locate molecu¬ 
lar clouds (see, e.g., Aravena et al. 2014) and reveal any possible 
molecular outflow component. 
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